Exterior optical cloaking and illusions by using active sources: a boundary element 

perspective 
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Recently, it was demonstrated that active sources can be used to cloak any objects that lie outside 
the cloaking devices [Phys. Rev. Lett. 103, 073901 (2009)]. Here, we propose that active sources can 
create illusion effects, so that an object outside the cloaking device can be made to look like another 
object. Invisibility is a special case in which the concealed object is transformed to a volume of 
air. From a boundary element perspective, we show that active sources can create a nearly "silent" 
domain which can conceal any objects inside and at the same time make the whole system look like 
an illusion of our choice outside a virtual boundary. The boundary element method gives the fields 
and field gradients (which can be related to monopoles and dipoles) on continuous curves which 
define the boundary of the active devices. Both the cloaking and illusion effects are confirmed by 
numerical simulations. 

PACS numbers: 41.20.Jb, 42.79.-e 



I. INTRODUCTION 



The classical wave scattering cross section of an ob- 
ject can be significantly larger or smaller than the ge- 
ometric cross section. 1,2 Recipes to achieve invisibility 
(zero cross section) are particularly intriguing, 3 nor- 
mally enabled by the concept of transformation media 
and artificial metamatcrials. 6,11-13 The correspondence 
between coordinate transformation and material parame- 
ters was noted nearly half a century ago 3,6 and such cor- 
respondence was explicitly formulated as the technique 
of "transform optics" to achieve invisibility by steering 
electromagnetic waves around a domain. 6,7 Similar ap- 
proaches to achieve invisibility were also proved mathe- 
matically for geometric optics 4,5 and in the quasi-static 
limit. 8 These invisibility devices typically work by steer- 
ing light around an object and the material shells need to 
encircle the object to be cloaked. 6 It was then proposed 
that "cloaking at a distance" can be achieved, 14 and the 
concept of cloaking can be extended to create arbitrary 
illusions. 15 However, these recipes are based on artificial 
metamaterials and usually have bandwidth limitations. 
Recently, cloaking by using active sources have been pro- 
posed, 16,18,19 which removes the requirement of metama- 
terials as well as the bandwidth limitation. Miller gave 
a detailed prescription to perform active source cloaking, 
and considered the case in which the sources encircle the 
cloaked domain. It was then shown by Vasquez, Milton, 
and D. Onofrei, 17 that exterior cloaking can be also re- 
alized by using several points (disks) of active multipolc 
sources placed around the object to be cloaked. This 
cloaking effect has been demonstrated numerically in a 
broadband fashion. 17,20 

In this paper, we will employ a boundary element for- 



mulation to show that arbitrary illusions can be achieve 
using simple active sources, i.e. sources of fields and field 
gradients (monopoles and dipoles), placed on continu- 
ous curves; and active source remote cloaking is a special 
case. When the active sources are tuned properly accord- 
ing to incoming waves, a nearly "silent" domain is created 
such that any objects can be hidden inside. At the same 
time, invisibility is achieved by reducing the "scattered" 
waves of the active sources to be almost zero on a vir- 
tual boundary which encloses the whole system. We also 
extend this scheme of external active cloaking to create 
arbitrary illusions, by tuning the active sources to pro- 
duce the same "scattered" waves on the virtual boundary 
as those scattered by the object chosen for illusion un- 
der the same incident wave. The physics of the cloaking 
and illusion effects by active sources can be understood 
clearly from the boundary element perspective. 21 



II. BOUNDARY INTEGRAL METHOD 

In this section, we demonstrate that the problem of 
two-dimensional (2D) active cloaking can be formulated 
by using a boundary integral equation 21 and be extended 
to create arbitrary illusions. A schematic figure of the 
configuration is shown in Fig. 1. The active sources are 
placed on boundary of the domains marked by the red 
lines (labeled as r o ). These active sources will generate 
fields so that any object inside a certain domain Q c will 
become invisible and the external observer (outside a vir- 
tual boundary Tf,) will see an illusion of another object 
inside. The active sources do not need to encircle the ob- 
ject to be cloaked, as shown in Fig. 1. For any incoming 
waves, these active sources generate opposite fields that 
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FIG. 1: (Color online) Schematic figure for the illusion devices 
containing exterior active sources on the boundaries F a (red 
lines) that create a cloaked region Q c which can hide arbitrary 
scatterers inside and at the same time make the total system 
response like another illusion scatterer Vi outside a virtual 
boundary IV When the illusion scatterer is just free space, 
invisibility is achieved. 



cancel the incoming waves inside the domain Q, c C R 2 to 
make the total fields inside J7 C almost zero. Any objects 
inside this "quiet" zone will be concealed as the objects 
experience no incident wave and thus no scattering will 
occur. 

At the same time, the active sources can be used to 
generate outgoing fields outside the boundary Tb that 
mimic those scattered from another object Vi under the 
illumination of the same incoming wave, rendering the 
whole system to appear like the object Vi for any ob- 
servers outside IV In other words, the wavefronts on 
r& of the scattering fields from the object Vi are recon- 
structed by the active sources. Thus, the active sources 
can create an illusion so that any object placed inside 
r c is transformed optically so that it looks like another 
object for observers outside IV If the fields generated by 
active sources cancel each other outside a virtual bound- 
ary (labeled as IV , then any observer outside Tb would 
not see the object inside T c as well as the active sources. 
In other words, the object inside J7 C is cloaked by the 
active sources. In this sense, invisibility is a special case 
in which the illusion object V, is just free space. 

The pre-designed object Vi, which is the illusion we 
want to create, can be changed on demand as long as 
we can compute the scattered fields of this object un- 
der the incoming probing waves. In this active source 
approach, the external illusion device does not need to 
be tailor made for the object and illusion, in contrast 
to the exterior cloak built with passive metamaterials, 15 
and there is no intrinsic bandwidth limitation. However, 
one needs to know in advance the probing wave, or one 
must set up sufficiently fast-responsive sensors to capture 
the information of the probing wave on the boundary of 

We will formulate the problem using the boundary el- 
ement approach. We note that the surface integral equa- 



tion (SIE) of the 2D Hclmholtz equation (V 2 + k 2 )4>(v) 
0, where k is the wave number, can be written as 21 

</>(r) = f ds[g(s,r)d a <j)(s) - 0(s)<9 n #(s, r)] 
r efi Jan 



(1) 



Here n is the outward unit normal vector on the bound- 
ary and <9 n represents the normal gradient. The 2D 
SIE tells us that inside a homogeneous domain 57, the 
scalar wave function </>(r) is completely determined by the 
fields and the normal derivatives on the boundary, con- 
nected by the Green's function ^(r, r') = \H^' (k\r— r'|), 

where Hq is the zeroth order of the first kind of Hankcl 
function. The counterpart of Eq. (1) for an open do- 
main which might also involves an incident field </> lnc (r), 
reads (f>(r)\ re ^ R 2_ n j = mc (r) + bc (r), where <f) bc is con- 
tributed from the following boundary integral 



^ bc (r) = - 



f ds[g(s,r)d n (t)(s) - (j>(s)d n g(s,r)] . (2) 
Jon 

Inspired by Eq. (2), one can construct additional fields 
(<f) hc ) in the domain R 2 — 57, by appropriately choosing 
active sources to create <p(s) and d n <p(s) on the boundary 
<957. In the Appendix, we show that such kinds of active 
fields can be generated by a variety of sources either lo- 
cated on the boundary or inside the boundary and these 
fields correspond to outgoing multipolc radiation fields 
from a perspective of an observer outside the cloaking 
devices. Using Eq. (2), one can determine each compo- 
nent of the multipole sources. Now, the issue is that 
whether it is possible to use such kind of active sources 
to construct fields which can cancel the incident field </> mc 
inside the "quiet" zone 57 c , and simultaneously mimic the 
scattered field in the region outside Tt (i.e., R 2 — 57;,), i.e. 



<^ bc (r) = 



-0 inc (r), 
<T(r), 
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(3) 



such that the total field is essentially zero inside f2 c , 
while outside 57;, the total field mimics the superimpo- 
sition of the incoming waves and the scattered wave of 
a pre-designed illusion object. As will be shown in Ap- 
pendix B, the perfectness of the cloaking and illusion 
effect depends on the number of active sources that we 
can afford to use. Employing a uniqueness theorem 22 , 
the conditions in Eq. (3) can be simplified. In general, 
fixing the values of </> on the boundary (Dirichlet bound- 
ary condition) can already guarantee a unique solution 
in the enclosed domain for the 2D Hclmholtz equation. 
Thus, the constraints in Eq. (3) can be replaced by 



bc (r) 



-0 inc (r), for rer c , 
<f c (r), for r E T b . 



(4) 



In addition to the forgoing conditions, there is a self- 
consistent condition of the boundary fields on r a , 21 



^(s) = <T c (s) - 



ds'{g(s',s)d n <f>(s') 
s')cVg(s',s)}, s,s'er o .(5) 
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in which, the integral is of a Cauchy principal value 
(CPV). This self-consistent condition comes from the 
continuity requirement when r approaches to the bound- 
ary from outside the cloaking device, i.e. 

Iim0(r) = ^(s), r ER 2 -n a ,seT a . (6) 

r— >s 

i 



The solution to these integral equations can be numer- 
ically determined using the boundary element method 
(BEM), 21 which is based on the SIE. BEM approximates 
the surface integrals by discretizing the surface T a = d£l a 
into N surface elements £ a on which the functions <^>(s) 
and ^( s ) = <9 n 0(s) are approximated as constants that 
represent the value of functions <j)(s) and t/j(s) across 
the entire element £ a , respectively. In other words, a 
local step function basis with regard to i a is used to 
expand </>(s) and ^(s) over the entire surface dfl, with 
the expansion coefficients denoted as <p a and ip a , where 
a = 1, 2, • ■ -N. From the perspective of BEM, these ex- 
pansion coefficients can be viewed as 2N active sources 
to be determined, and Eq. (7) represents the conditions 
to determine these 2N unknowns. Equation (7a) gives 
N constraints. Besides these, one can choose N c sample 
points on T c and iV& points on Tb and then get a total of 
N c + Nj, + N constraints to determine the 27V degree of 
freedom <j> a and ijj a ■ We can then establish the following 
linear equations 
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where H(, a , G& a and H ca , G ca represent the interacting 
matrices relating the "source points" on r a to the field 
points on Tfc and r c , and have elements defined as 

H a/3 = / d n g{s,r a )ds , (9a) 
J tp 

G a f3 = g(s,r a )ds (9b) 

Jtp 

r a 6 r c ur 6 ,a = 1,2,... ,N C + N b , 
If, cT a ,/3 = l,2,...,N. 



In summary, the active fields are determined by the fol- 
lowing integral equations, 



= <T c (r), rer a , (7a) 



</> sc (r), rer b , (7b) 



-<T c (r), r G r c . (7c) 
I 

whereas, H™* and G™* represent the self-consistent con- 
ditions imposed on <fi and d n (f>, with elements defined as, 

= - [ d n g(s,r a )ds, (10a) 

G' Q n | = / g(s,r a )ds (10b) 

Jig 

e a ,lpCT a ,a,0 = l,2,...,N, 
r„ is the center of £ a . 

On the right hand side of Eq. (8) , 3>J, nc denotes the incom- 
ing probing wave fields lnc sampled at the inner quite 
zone boundary T c = dfl c , <b s b c denotes the sampled scat- 
tered fields on the outer boundary Tf, = dflb that would 
have been scattered by the object Vi under the illumi- 
nation of the same probing wave </> lnc . $ 1 a Ilc denotes the 
incoming wave fields on the boundary of the cloaking 
devices. Thus, the fields outside Tt (the dashed-line in 
Fig. 1) approach those of 4> lnc scattered by the object Vi 
as the number of sampling points(N) increases and the 
discrepancy decreases if N increases. If we set €>^ c = 0, 
we achieve the active external invisible cloaking. Cloak- 
ing is thus a special case of illusion in this formulation. 
<fr a and ty a are both ./V-dimensional vectors represent- 
ing the total field and field gradient on T a , which require 
active sources to generate. For simplicity, in our numer- 
ical calculations, we set N c + Nb + N = 2N, so that 
the matrix in Eq (8) is a square matrix. The linear sys- 
tem of equations is solved using the LAPACK subroutine 
ZGESV. We can see that BEM offers a physically trans- 
parent way of deriving the active sources needed to do 
remote cloaking and illusion and offers a straight forward 
numerical recipe in calculating those sources. 

Our forgoing discussions are restricted to non-radiating 
objects. If the object inside f2 c itself is a radiating source, 



\(j>{v) + / {g(s, r)d„0(s) - cj)(s)d n g(s, r)} ds 
z Jon a 



on a 



{g(s, r)<9 n c/>(s) - (j)(s)d n g(s, r)} ds = 



{g(s, r)d a <f>(s) - (/)(s)d n g(s, r)} ds 
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we should add extra active sources to cancel the radiated by 
field radlatc outside IV This extra term is determined 



J 



Z Jdfla 



Radiate (r)) reT ^ 



(11a) 



/ {5(s,r)a n cxtra (s)-0 cxtla (s)9 n g(s,r)}ds = -0 radiatc (r), r e T 6 , 
"'SO, 



(lib) 



/ {.g(s,r)a n cxtla (s)-0 cxtl ' a (s)a n g(s,r)}ds = 0, r e T c 



(lie) 



The Green's function g might be different since the radi- 
ating field may be of a different frequency <D. Then the 
total active field should be 

4> a (r, t) = 4>(r) exp [-iut] + cxtra (r) cxp [-iwt] , (12) 

where 0(r) is the solution of Eq. (7), and GXtra is the so- 
lution of Eq. (11). We note that the condition in Eq.(llc) 
is necessary since the radiating object may also be a pas- 
sive scatterer. 



III. NUMERICAL SIMULATIONS 

A. Remote Cloaking Effect 

Here, we show numerically that external cloaking is 
achieved by simply setting & s b c = in Eq. (8). The choice 
of plane wave exp(ik-r) as the incoming source is just for 
simplicity but the formulation works for other forms of 
incident wave. The configurations of the cloaking devices 
are shown in Fig. 2 for two kinds of source arrangement. 
In the left panels, active sources are placed on three cir- 
cles arranged as shown in Fig. 2(a) and Fig. 2(c). In the 
right panels, we show a case in which the active sources 
are placed on two crescent-shaped curves [Fig. 2(b) and 
Fig. 2(d)]. In both cases, we have chosen N = 900, 
N c = 300, and Nf, = 600. Employing the scheme de- 
scribed in the preceding section, we can achieve an ap- 
proximate solution numerically. We see from Fig. 2 that 
the field inside the quiet zone is essentially zero and here 
is no scattering. The numerical solutions only ensure 
that Eq. (7) is correct in a finite number of points. To 
quantify the overall quality of the solution, we consider 
the following error functions (measured with L 2 norm) r e _g 
defined on the two circles Tf, (r=20) and T c (r=2) and 
and inside the quiet zone f2 c , 



Err(r c ) 



Err(r 6 ) = 



& n >(s)-^ c (8)| 2 d5 



c (s)| 2 ds 



fan c 



(13b) 



(13c) 



JdSl b 



= (s)| 2 ds 



(13a) 



Each integrand is numerically evaluated at a set of 40,000 
points sampled in the corresponding integration domain. 
For fixed circle sizes of and r c , the errors depend on 
the choice of N, Nb and N c , as well as the frequency of 
the incoming wave. A detailed discussion on the depen- 
dence of these errors on the parameters can be found in 
Appendix B. It is seen that the error decreases as we 
increase N. In other words, we are able to achieve bet- 
ter cloaking effects if we can control the boundary fields 
more precisely. The field patterns are presented in Fig. 2. 
Figs. 2(a) and 2(b) show the total fields which are the 
superposition of the incoming plane wave and the ac- 
tive fields generated by the "active sources" [solutions to 
Eq. (8)] placed on the white dotted-lines. The total fields 
</>(r) outside Tb (marked by the black dashed- line), a cir- 
cle of r = 20 units, are almost the same as the incoming 
plane wave, with a discrepancy Err(Fb) = 4.62 x 10~ 13 
in Figs. 2(a) and 2(c). At the same time, we achieve 
a "quiet" zone f2 c (bounded by r c as marked by the 
green lines) within a circle of r = 2 inside which the 
total fields almost vanish, with Err(r c ) = 1.14 x 10~ 12 
and Err(fi c ) = 1.03 x 10~ 12 . Figs. 2(c) and 2(d) show the 
corresponding "scattered" fields </> sc (r) = 0(r) — ^> lnc (r), 
which are exactly the fields created by the active devices. 
Concomitantly, </> sc (r) vanishes outside Tf, and is the re- 
verse of </> lnc (r) inside the "quiet" zone f2 c . The strength 
of the fields on the boundary T a is in the order of 100 
e.g., - 160 in Figs. 2(a) and (c), ~ 80 in Fig. 2(b) and 
(d)], and can be easily achieved physically. Vasquez and 
co-workers 17 proposed that three disjoint circular disks 
are needed to perform remote active cloaking. Here, we 
see that from Figs. 2(b) and 2(d) that a non-circular 
cloaking device comprising two simply connected regions 
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FIG. 2: (Color online) Example of active source external cloak. (a),(b) show the total fields, and (c),(d) show the "scattered" 
fields. The sources are arranged on the boundaries (white solid line) of three circles in the left panels while are in two 
crescent-shaped curves in the right panels. Here, the incoming plane wave is of wavelength A = 3.0. 



can also achieve the similar cloaking effect. It is further 
shown (figure not presented here) that active sources on 
one simply connected cloaking device can also achieve a 
high degree of invisibility. 



B. Illusion Effect 

Next, we demonstrate an illusion effect such that what- 
ever objects placed inside the "quiet" zone fl c will appear 
like another object, which is chosen here to be a banana- 
shaped dielectric object with refractive index n = 2.32. 
Figs. 3(c) and 3(f) show the total and scattered fields of 
such a "banana" under the illumination of a plane wave 
cxp(ik • r). We set the control boundary (marked 
as the dashed-line) as a circle with radius r = 20 out- 



side which the illusion shall be observed. This requires 
one to set €>^ c in Eq. (8) as the scattered fields due to 
the pre-designed "banana" [see Fig. 3(c)] onTj,. A total 
of N = 900 sample points (with N b = 600, N c = 300) 
are used here in the numerical calculation, and $J, nc in 
Eq. (8) is the function exp(ik • r) sampled over a circle 
of r = 5 (the boundary of the "quiet" zone fl c as marked 
by the green solid circle). In Figs. 3(b) and 3(e), we 
conceal an apple-shaped dielectric object (n = 4) inside 
the "quiet" zone. Figs. 3(a) and 3(d) show the total and 
scattered fields of the "apple" under the illumination of 
the plane wave. We see in Figs. 3(b) and 3(e) that, after 
turning on the active sources, the total system responses 
to the incoming plane wave in a way as if a "banana" is 
placed in the cloaked region [Figs. 3(c) and 3(f)]. The 
field discrepancy on the circle of r = 20 is 1.17 x 10 ~ 6 , 
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measured with L 2 norm as, 
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Err iH (r 6 ) 



^ o j0(s)-0 i " c (s)-^(s)| 2 d S 
§ an W™{s)\Hs 



, (14) 



where <fi sc is the scattered field of the pre-designcd "ba- 
nana" under the illumination of the incoming wave. 
Fig. 3(e) shows that the active sources construct a neg- 
ative counterpart of the incoming plane wave to create 
the "quiet" zone f2 c , so that the total field surrounding 
the "apple" is almost zero. We note that we can put any 
non-radiating object inside O c without affecting the to- 
tal scattering pattern outside T b [see Figs. 4(b) and 4(c)], 
since all the fields inside f2 c are almost zero as is shown 
in Fig. 4(a) [The discrepancies: Eri J "(r c ) = 6.15 x 1CT 7 , 
Err in (f2 c ) = 1.11 x 10 -6 ]. In fact, any passive object in- 
side fl c does not "talk" to the other parts of the world, 
rendering this illusion device workable for multiple and 
arbitrary objects. In principle, by changing the active 
sources, we can let observers outside T b see whatever we 
want them to see. 
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Appendix A: Possible ways to build up the active 
sources 

The boundary element method gives the field and field 
gradients on the boundary of the active cloaking device. 
It does not give directly the sources that give rise to those 
fields. There are obviously many different ways of arrang- 
ing sources within the boundary of the active devices to 
generate the necessary fields. In this Appendix, we give 
a few examples. 



Using monople sources and dipolar sources 
located on the boundaries 



C. Cloaking a radiating object 

Here, we demonstrate the cloaking effect of a radiat- 
ing object. For simplicity, we assume that ^ radlato = 
10-h[ (kr) cos 9 exp[— iurt] and there is no incoming wave 
from the outside. The simulation results are shown in 
Fig. 5, with k = it. Here, we use the same parame- 
ters N, N b , N c as those in section III A and section IIIB. 
From Fig. 5(b), we see that the total field pattern near 
the dipole remain the same with Fig. 5(a). The cloak- 
ing device creates active fields which cancel the radiating 
fields outside r& [see Figs. 5(b) and 5(c)] and do not af- 
fect the field pattern inside fl c . The active source in 
this example is just the c/> extra (s) which is the solution of 
Eq. (11). Considering the linear superposition property, 
one can add together the solution of Eq. (7) and ^ extra 
to achieve cloaking or illusion effects for this specific ra- 
diating object. 



IV. CONCLUSION 

In conclusion, we applied a boundary element method 
to demonstrate both external cloaking and illusion ef- 
fects using active sources on continuous curves for the 
2D Helmholtz equation. The scalar wave formalism ap- 
plies to both acoustic waves and electromagnetic waves 
in two dimensions. This approach works for arbitrary 
objects and there is no intrinsic bandwidth limitation. 
The limitation of this type of remote active cloaking and 
illusion is that it requires the prior knowledge of the in- 
coming wave or the availability of sensors that can detect 
the fields quickly enough on the boundaries and of active 
sources that can respond fast enough. 16 



We note that the active fields are given by, 

activc (r) = _ / ds [ giS:r)dMs) _ a ( S )d nff (s,r)] , 
JdQ a 

(Al) 

where </> a and d n (j) a are the solutions of Eq. (7). The first 
term represents a set of monopole point sources, with the 
corresponding strengths as follows 



Pi = -9 n 4> a (si)Sai , 



(A2) 



where is the center of the i-th boundary element, and 
5<Xi is the length of the corresponding element. The sec- 
ond term represents a set of dipole sources on the bound- 
ary, since 

d n g(s,r) = ^W^r-sD^p 
= !*!_!#)(*[, _ B | )e " 

+^^H[ 1 \k\r^s\)e^, (A3) 

which are exactly dipole sources located at Sj. Here, 6' is 
the direction of the unit normal n and 9 is the direction 
of r — s. We have used the recursion relation of Hankcl 

functions = -H^ix). 

The corresponding strengths of the two types of dipole 
sources are 



„(*), forff^fclr-SiDe-tA*) 



It can be easily verified that this set of point sources and 
dipole sources can exactly generate the necessary active 
fields [see the schematic configuration on Fig. A 1]. 
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FIG. 3: (Color online) Optical illusion effect by active sources, (a), (b) and (c) show the total fields, (d), (e) and (f) show the 
"scattered" fields (0 tot — <^ mc ). (a) and (d) are for an apple (e = 16.0, fx = 1.0). (b) and (e) are for an illusion device with an 
apple-shaped object [identical to (a) and (d)] concealed inside Q c . The active sources are placed on the three circles (white 
solid line), (c) and (f) are for a banana-shaped dielectric object (e = 5.0, fi — 1.0). The active sources give the same scattering 
pattern as the "banana" outside the black dotted curve. 



2. Using double layers of point sources 

We can also use only point sources to build up the 
same necessary active fields 16 . Instead of positioning the 
sources on the boundary T a , one can place point sources 
on the outer boundary S out and the inner boundary Si„ 
[see the schematic configuration on Fig. 7] , i.e. we place 
point sources at the following points 

s ?«* = Si + | ni , s f = s ,-| ni , i = l,2,..,iV(A5) 

where is the center of the i-th element and is the unit 
normal of the i-th element. The corresponding strengths 
of the point sources are, 



It can easily seen that the field generated by the point 
sources above is 

JV 

0(r) = J2[p( s ° Ut )9(s° Ut ,r)+P(sT)9(sr,r)_ 

i 




ds[g{s,r)d n (j) a (s) - <j) a (s)d n g(s,r)] 



+0(h 2 ) . (A7) 

As long as the distance h between the outer and the inner 
boundary is sufficiently small, one can build up the active 
field as close as we like to that of Eq. (Al). 

3. Using multipole sources placed at the centers of 
the cloaking or illusion devices 



P(s° Ut ) 



~h . 

5 a 
h 



(A6) 



Instead of positioning monopole or dipole sources di- 
rectly on the boundaries of the cloaking or illusion de- 
vices, one can put multipole active sources inside the 
cloaking devices Q a to construct the required boundary 
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FIG. 4: (Color online) Optical illusion effect by active sources, (a), (d) show the total field pattern of an illusion device, (b) 
and (c) are patterns of the total fields for the same illusion device but with different objects concealed inside fl c [(b) is a triangle 
e — 2.25, (i = — f .0, (c) is a PEC slab], (e) and (f) are patterns of the total fields for the objects identical to those concealed 
inside the "quiet" zones in (b) and (c). Any object can be concealed inside the "quiet" zone where the total fields are almost 
zero. 



fields (j) and <9 n 0, and then control the fields outside tt a . 
For a circular domain, it is convenient to place the multi- 
pole source at the center. In the following, wc will discuss 
a way to determine the corresponding strength of each 
component from the boundary fields obtained by BEM. 
We note that the active fields are given by, 

^activc (r) = _ I ds [g {Sir)dMs) _ a ( S )a n5 (s,r)] . 

Jon a 

(A8) 

For the domain outside the cloaking device, | s — r c | < 
r — r c | (r c is the center of the cloaking device), we can 
expand the Green's function as, 

oo 

g(r,s) = H^(k\r-v c \)J m (k\s-v c \)exp[im(A9)}, 

m— — oo 

|s-r c | < |r-r c |, (A9) 

where i?™' is the m-th order Hankcl function of the first 
type and AO is the angle between r— r c and s — r c . Substi- 
tuting this into Eq. (A8), one can get the corresponding 



multipole sources. For example, for the case shown in 
Fig. 2(a), one can use multipole sources locating at the 
centers of the three circles of the circles r Ci , i = 1, 2, 3 to 
generate the required fields and field gradients, i.e, 

3 oo 

activc (r) =J2 a ^ H rn ^ ~ ^ I) ex P [zmfl] , 

2=1 rn— — oo 

(A10) 

where 9 is the direction of r — r Ci , and a m .i is the strength 

of the multipole source locating at r Ci , which can 
be obtained by the following integrals, 

a m ,i = f dscxp[im9']l d n J m (k\s - r Ci |) cj) a (s) 
Jan a . L 

-J m (fc|s-r Ci |) d„0 a (s)}, (All) 

where i = 1,2, 3, J m (x) is the m-th order Bessel function 
and 9' is the direction of s — r Ci . 

The multipole expansion of the Green's function 
demonstrates that the two kinds of boundary fields and 
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FIG. 5: (Color online) Cloaking a radiating object. (a)The radiating field pattern of 10 • H^^kr) cos 8 exp[— iuit], where k = it. 
(b) The total field pattern with the sources on the boundary of the cloaking device actives(three circles) radiating according 
to the solutions of Eq. (11). (c) The field generated by the cloaking device. 





Pout 



Pin 



FIG. 6: Generating active field by using monople sources and FIG. 7: Generating active field by using point sources located 
dipolar sources located on the boundaries of the cloaking or on double layer boundaries, 
illusion devices. 



their gradients (4>(s) and d n cf)(s)) contribute the same 
form of multipole radiation fields H m {kr) exp(im0), with 
their strengths determined by the surface integral. In this 
sense, 4>(s) and d n <fi(s) are on the same footing in generat- 
ing the active field in the open domain(i.e. R 2 — fl a ). This 
indicates that, in Eq. (A8), there are multiple choices of 
the boundary fields cf> a and d n (j> a , for the same <^ actlve (r) . 
This is reasonable, as one can easily verify that, 

/ d S [0 fc (s)a n5 (s,r)- 5 (s,r)a n( /) fe (s)] =0, 

JdQ a 

for r e R 2 - n a (A12) 

if 4> k (r) satisfies the Helmholtz equation 

(V 2 + k 2 )cf> k (r) = 0, r e Q a . (A13) 

If a {(f> a (s) , d n (f> a (s)} pair can give the required active 
sources, any solution of the form {0 a (s)+0 fc (s), <9 n a (s) + 



d n <j> k (s)} gives the same active fields. Then a question 
arises: Is the solution obtained by Eq. (8) unique? In 
fact, we have imposed the continuity boundary condition 
across T a , 

lim^r) = 0( s )| seaOa . (A14) 

r— 7-s 

This automatically removes the extra degree of freedom, 
since any extra term cf> k (s) makes the field discontinuous 
when taking the limit to the boundary. Physically, if we 
are using multipole sources located inside the cloaking or 
illusion devices, the continuity boundary conditions must 
be satisfied. 

However, as we know from Eq. (A12), the extra term 
4> k does not influence the active fields outside the cloaking 
device. Therefore, instead of imposing the boundary con- 
tinuity conditions, one can first find solutions of the form 
{^ Q ,0} or {0,d n <p a }, since both <j> and d n 4> can indepen- 
dently provide the required active sources "mathemati- 
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cally". In fact, physically, {4> a , 0} corresponds to placing 
dipoles on the boundaries and {0,d n <p a } corresponds to 
placing monopoles on the boundaries [see Eqs. (A2) and 
(A4)]. It is seen that one can use only the monopoles 
or only the dipoles placed on the boundaries to obtain 
the required fields. For the {4> a ,0} type of solution, the 
active fields can be expressed as, 



r cUve {r) = - / g{s, r)d n a (s)d S , r G R — f2 
Jon a 

(A15) 

Under this circumstance, the conditions for cloaking that 
the active fields should satisfy are, 



on a 



g(s,r)d n 4> a {s)ds 



3 (r), for r e T c , 
0, for r e T b . 

(A16) 



One can also find the corresponding multipole sources, 
with the following strength coefficients 

a m .i = ~f ds cxp[im0'] J m (k\s - r c J) <9 n 0(s) . 
Jon ai 

(A17) 

The set of a m ,i should be the same with that obtained 
from the first method in Eq. (All). Please note that, 
in this case, we does not impose the boundary continuity 
conditions, therefore, if one takes the following limit from 
outside the cloaking device, 

<^ s -(s)| se0 n = lim0(r), (A18) 

r— >s 

d n r hys -(s)\ sedn = lim3 n 0(r), (A19) 



same as the incoming wave outside the circle of r = 20.0 
and nearly zero inside the quiet zone inside the circle of 
r = 2.0, which confirms that for all the practical pur- 
poses, objects concealed inside the "quiet" zone are un- 
detectable outside of IV To get a better picture of the 
performance of the cloaking as TV increases, we provide 
some representative field patterns of the three schemes 
in Fig. 9. As can be clearly seen, when TV is not large 
enough, there will be a lot of evident "tails" reaching out 
from the boundary of the cloaking device and the field 
values are large, while for increasing TV, the "tails" be- 
come less conspicuous. Meanwhile, the field patterns of 
the three different schemes look indistinguishable with 
each other as TV increases. 



generally, ^> phys -(s) is different with <p a (s). However, in 
principle, the fields {</> phys -, d„<T> phys } should be the same 
with the {4> a , d n 4>a] acquired by imposing the boundary 
continuity conditions if the number of sample points N is 
large. This method can reduce the dimension of the ma- 
trix by half. We present the numerical simulation results 
of different formalisms in Appendix B. 



Appendix B: Convergence discussion 

The numerical solution of Eq.(8) only ensure that 
Eq. (7) is satisfied for a finite number of points. Here, 
we consider the convergence of the solution as the num- 
ber of sampling point increases by calculating the errors 
on the two circles r c and Tb. The errors are measured 
as Eq. (13) in the text. As shown in Fig. 8, the error 
decreases very quickly as TV increases (note that the ver- 
tical axis is in log scale). When TV is greater than 300 
(this number, however, depends on the wavelength and 
the geometry of the cloaking device), the total error is of 
~ 10~ 14 . This indicates that the field are essentially the 
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150 200 

N 



150 200 

N 



301) 



(a) r c , r = 2.0 



(b) F b , r=20 



FIG. 8: (Color online) Errors on (a) r c : r — 2.0 and Ft: r — 20.0, The wavelength A = 10.0. The errors are measured with the 
L 2 norm, defined in Eq.(13). 
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